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Jorge Castineiras 

Institute de Fisica Teorica, Universidade Estadual Paulista, Rua Pamplona 145, 01405-900, Sao Paulo, Sao Paulo, 

Brazil 

Luis C. B. Crispino 

Departamento de Fisica, Universidade Federal do Para, Campus Universitario do Guama, 66075-900, Belem, Para, 

Brazil 

George E. A. Matsas 

Institute de Fisica Teorica, Universidade Estadual Paulista, Rua Pamplona 145, 01405-900, Sao Paulo, Sao Paulo, 

Brazil 

Daniel A. T. Vanzella 

Physics Department, University of Wisconsin-Milwaukee, 1900 E Kenwood Blvd., Milwaukee, WI 53211 

We analyze free elementary particles with rest mass m and total energy E < mc? in the Rindler 
wedge, outside Reissner-Nordstrom black holes and in the spacetime of relativistic (and non- 
relativistic) stars, and use Unruh-DeWitt-like detectors to calculate the associated particle detection 
rate in each case. The (mean) particle position is identified with the spatial average of the excitation 
probability of the detectors, which are supposed to cover the whole space. Our results are shown to 
be in harmony with General Relativity classical predictions. Eventually we reconcile our conclusions 
with Earth-based experiments which are in good agreement with E > mc? . 
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I. INTRODUCTION 

The standard theory of quantum fields uses the fact 
that the Minkowski spacetime is maximally symmetric. 
The linear three momenta (fc^,A:^,fc^) associated with 
the translational isometrics on the spacelike hypersur- 
faces t — const constitutes a suitable set of quantum 
numbers to label free particles, where we are assum- 
ing here that it, x, y, z) are the usual Minkowski co- 
ordinates. In this simple case, the dispersion relation 
E = huj = -y/jkcp + c"* imposes a simple constraint 
between the particle mass to, momentum k and energy E, 
and, thus, free particles with well defined linear momenta 
must have total energy E > mc^. Because Earth-based 
particle experiments assume in general the detection of 
asymptotic free states in Minkowski spacetime, the pos- 
sibility of measuring particles with total energy 



E < mc 



(1.1) 



is usually neglected. Moreover, in the classical context of 
General Relativity, the detection in loco of point particles 
satisfying Eq. (LI) by direct capture is ruled out by the 
fact that an observer with four-velocity intercepting a 



particle with four-momentum — mv^ ascribes to the 
particle an energy E — tow^m^ > mc^. On the other 
hand, it is well known that, unlike standard Quantum 
Field Theory in flat spacetime, the field quantization car- 
ried over curved (stationary) backgrounds does not lead 
in general to any dispersion relation for the energy and 
other quantum numbers, avoiding thus the flat spacetime 
constraint E > mc^. This can be understood by recall- 
ing that the concept of point particle finds no room in 
the context of Quantum Field Theory. Two main ques- 
tions are raised, then, by the discussion above: Given 
a stationary spacetime (i) what is the minimum energy 
i^min = hujmin allowed for a particle ? and (ii) what is the 
probability density associated with the detection of parti- 
cles with E e [i^niin, "^c^) at different space points? 

We analyze the questions above in the context of Quan- 
tum Theory of Linear Fields in Curved Spacetimes where 
the normal modes associated with our particles are seen 
as free. In order to define clearly what we mean by "ob- 
serving a particle" we use Unruh-DeWitt detectors (en- 
dowed with an internal structure defined by a density of 
states). Our results are shown to be in agreement with 
General Relativity predictions concerning the position of 
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particles with E < mc^ as defined by the same fiducial 
observers who proceed to the field quantization and with 
Earth-based experiments which usually assume E > mc^ . 

The paper is organized as follows. In Section || we con- 
sider a two-dimensional Schwarzschild black hole which 
maintains a close relationship with the Rindler wedge 
where a full analytical investigation can be carried out. 
This section is important as a "theoretical laboratory" 
for the next ones and to provide a better understanding 
of some phenomena as, e.g., the decay of miiformly ac- 
celerated protons [Q . This is so because the proton inter- 
action with massive particles possessing energy E < mc^ 
plays an important role according to coaccelerated ob- 



servers with the proton. In Section [II we investigate the 



questions (i) and (ii) outside Reissner-Nordstrom black 
holes. We discuss, in particular, the agreement of our 
results with General Relativity. In Section IV we con- 



sider relativistic (and non-relativistic) stars. This case is 
qualitatively different from the earlier ones because the 
energy spectrum of particles with E < mc^ is discrete. 
This fact is a reflection of the nonexistence of event hori- 
zons in the star spacetime. We show eventually that for 
Earth-based experiments i^min ~ mc? , as expected. Our 
final remarks are made in Section ^ We will assume 
hereafter natural units, fi = c = G = 1, unless stated 
otherwise. 



II. DETECTION OF MASSIVE PARTICLES IN A 
TWO-DIMENSIONAL BLACK HOLE-LIKE 
SPACETIME 

Let us begin considering the line element of a two- 
dimensional Schwarzschild spacetime: 



(1 - 2M/r) dt^ - (1 - 2Af/r)"^ dr"^ 



(2.1) 



Eq. (2.1) can be seen as describing a two-dimensional 
black hole with mass M . Close to the horizon, r « 2M, 



Eq. {2A_) can be written as 



ds^ = [p/AMfdt^ - dp^ 



(2.2) 



where p{r) = yJ^M{r — 2M). (Note that in these coor- 
dinates the horizon is at p = 0.) Line element (^^) is 
associated with the Rindler wedge (which is a globally 
hyperbolic spacetime) provided that < p < +oo and 
— oo < t < +00. The advantage of considering space- 
time (2.2) rather than ( ^.l| ) is three-fold: first it has 
the main relevant features (for our purposes) of the two- 
dimensional Schwarzschild spacetime (the fact that they 
differ asymptotically will not be important at this point). 




FIG. 1. We plot the detection rate Tuiq as a function of 
the uio/m ratio for observers at different points po- For ob- 
servers far away from the horizon, there is a strong detection 
damping. 



second it provides a better understanding of some phe- 
nomena occurring in uniformly accelerated frames and 
third it allows a completely analytic discussion. 

Let us now consider in this spacetime a free scalar field 
^(x) with mass m. The positive-frequency solutions of 
the Klein-Gordon equation (□ + m'^)uui{x'^) = can be 
written as u^{x^^) = V'w(p)e^*"* where 



dx"^ 



+Krv'.(p) = ^v.(p) 



(2.3) 



Here x = 4Mln(p/4M) and V;^^ = (4M)~ 'p- 
that the effective potential grows unboundedly at the in- 
finity. By solving Eq. (|2.3|), we obtain 



u^{x^') = V(4M/7r2) sinh(47rAfw) KuMu:{mp)e- 



(2.4) 



where K^, (x) is the modified Bessel function and we have 
orthonormalized u^{x) according to the Klein-Gordon in- 
ner product Q- [Q. Note that w € (0, -l-oo), i.e., there are 
massive Rindler particles with arbitrarily small energies. 

Next, the field is expanded in terms of positive- and 
negative-frequency modes as usually: 



$(a;) = / duj [d^u^{x) + H.c] , (2.5) 
Jo 

where the annihilation and creation operators satisfy 



5{uj — Lj'). Since dtu^^ 



-iuju^, the fidu- 



cial observer with respect to whom the quantization is 
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performed is the one at p = po — 4M, whose proper time 
is t [see Eq. (2.2)]. The Rindlcr vacuum |0) is defined 
from auj\0) — 0. 

Now, we introduce an Unruh-DeWitt detector |^ de- 
scribed by a locahzed monopole m{s) with proper time 
s and worldhne = z^(s). Let H be the detector free 
Hamihonian acting a,s H\E) — E\E) on the detector en- 
ergy eigenstates \E) and m{s) — e'^*?fi(0)e~'^''. We 
wiU denote by \Ea) the detector ground state and as- 
sume Eg = 0. The excited states of the detector wih 
be ruled by some (normalized) density of states Peq (E) 
peaked at = E'o and satisfying 



-t-CXD 



dEPEoiE) = 1 



(2.6) 



The simplest choice is /3eo{E) = 6{E — Eq) which char- 
acterizes a single excited-state detector. This will be 
enough in the black hole cases but not in the relativis- 
tic star one where we have to deal with particles with 
discrete energy spectra. For the later case it is more con- 
venient to consider 



/3e„{E) = {n/Eo) [e {E-Eo + Eo/2n) 
- e{E-En- En/2n)] , 



(2.7) 



where n = const. and Q{x) is the step function. 

We note that Eq. (^) satisfies Eq. (2^) and the useful 
property j3cEo{cE) — c~^Peq{E) for c € R. Moreover, 
we can recover the single excited-state detector in the 
n +00 limit, since lim„^+oo /3_B(, (i?) ~ S{E — Eq). 
Once the detector is defined, we couple it to a massive 
scalar field $(x'^) through the interaction action 

/ + 00 
dscom{s)^[x''{s)], (2.8) 
-oc 

where Cq is a small coupling constant. 

Let us now ask to our fiducial observer what is the 
total probability per (detector) proper time Ti^{pci) = 
Puj{pd) / s^d^ of detecting a particle at some point pd with 
energy lo. The excitation amplitude (at the tree level) 
At°^ = (0| (g) {E\Si\Eg) «) \uj) associated with the particle 
detection can be shown to be 



At''' = 4 Co ^M^hM^Mj) Ki.MUmpd) 

X S{E - UlLj/pd) , (2.9) 

where the detector selectivity was chosen such that 
{E\m{0)\EG) = 1. The detection rate is thus given by 

(3eo{E) / dcj'lAt'i'l'FU^') 



2 r-l-oo 

T^{pd) = — dEpEoiE) 



2c, 



— ^ smh{T:E(iPd)pdKfE„p^{'rnpd)F, 



Eopd 
AM 



(2.10) 



where we have chosen the density of states f3Eo{E) = 
6{E — Eq) and Fi^{uj') characterizes a mixed state peaked 
at the particle energy uj with the property that Ft^(uj) — 
no — const, for every lo. (One can avoid the introduc- 
tion of mixed states by considering the detection of wave 
packets, as shown below.) Now, we carefully adjust the 
detector at each point pd to maximize the detection prob- 
ability ( 2.10| ) . This is achieved by properly tuning its en- 
ergy excitation gap: Eq — AMco/pd (note that Eq and lo 
are related by a red-shift factor, as it should be). Hence 
Eq. ( |2.10D becomes 



2c^ 

Tu^ipd) = —no smh{'i7rMLo)pdKlMJrnpd) ■ (2.11) 

TT 

Now, we may wonder what is the probability rate T^^g 
that a massive state jwo) be measured by our observer at 
Pq in the particular case where the detector is carried by 
our experimentalist, i.e., pd = po (and we recall that in 
the Rindler wedge M = po/4). In Fig. (^ we plot 



—no smh{TrpoUjQ)pQKf 

IT 



■ipaLOo 



(mpo) , (2.12) 



as a function of LOo/m. We clearly note that the farther 
away the observer is from the horizon, the steeper the 
detection damping for coo/m < 1. 

Next, let us define from Eq. (2.11) the normalized 
probability density 



dV^/dpd = T^{pd)/ / r^{pd)dp'd 



-\-oo 



(2.13) 



{dVu:/dpd)dpd is the probability that a particle with en- 
ergy LO be found between pd and pd + dpd. We see from 
Fig. (|^) that observers far away from the horizon will only 
be able to interact with the "tail" of the "wave functions" 
associated with particles with small co/ni. 



The physical content carried out by Eq. (2.13) can 
be reproduced by considering wave packets (rather than 
mixed states) as follows. Let us represent a particle with 
typical energy to through a wave packet defined as 



-f oo 



duj'GU^') al\0) 



(2.14) 



where G^^^lo') is a peaked fmiction at lo' = lo and 
J^°° duj'\G^{Lo')\'^ = 1 in order that = 1. The 

total probability for detecting \(f>^) (at the tree level) 
at some pd between the constant-coordinate-time hyper- 
surfaces t — const — oo and t — const — > 4-oo is 
^r(P<i) = m® {Eo\Si\Eg) ® 10^) P , where we should 
tune the detector as before, Eq = AMuo/ pd^ to maximize 
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FIG. 2. We plot the probability density dVui/dpd for differ- 
ent u)/m ratios, where we have assumed Mm =1/4. We note 
that the smaller the uj/m ratio, the closer to the horizon (in 
average) the particle lies, where the "gravitational potential" 
is stronger. 



FIG. 3. {pd) is shown to be smaller than p™''" = AMu/m in 
the "high-frequency" regime u > (4M)~^ (i.e., at the right of 
the vertical broken line) just as expected. (We have assumed 
here mM = 10 but this nice agreement is verified for any 
mM.) 



its detection probability. In order to obtain the prob- 
ability between the constant-proper-time hypersurfaces 
s = const — > —oo and s = const — > -l-oo, we must multi- 
ply both sides by the red-shift factor AM/ pd'. 

^^(Pd) = 4c2|G^(a;)|2 sinh(4^Afa;)p<i KImu. {mpd) , 

(2.15) 



where Guj{uj) — const. Note that Eqs. ( ^.15 ) and (2.11) 
differ only by a (dimensional) const ant f actor. Thus by 
replacing T^{pd) by V^^{pd) in Eq. (|2.13|) , we obtain the 

same probability density dVu/dpd- 

Now, in order to interpret Eq. ( |2.13|) on General Rela- 
tivity bases, let us first consider a row of detectors each 
of them lying at different pd and define the average de- 
tection position 



dpd Pd dVuj/dpd ■ 



(2.16) 



By using Eq. (2.13), we obtain 

7rtanh(47rMu;)(64M2w2 



iPd) 



1) 



64mMw 
TT Mlo /m (w S> a) 



(2.17) 



where a = 1/ 4M is the proper acceleration of the fiducial 
observer. 

Now, from General Relativity, a classical particle with 
mass m lying at rest at some point pp has, according 



to our fiducial observer at po — 4Af, total energy uj — 
mpp/AM. By considering that the particle may have in 
addition some kinetic energy, the total energy would be 
to > mpp/AM. By inverting this equation, we obtain 



Pp 



< 



AMuj/m , 



(2.18) 



d] < Pp , 



which is expected to agree with (pd), i.e., {p, 
at least in the "high-frequency" regime to a (where 
the quantum and classical behaviors may be compared) . 
This conclusion is indeed in agreement with Eqs. (2.17) 
and (2.18) (see also Fig. ^). The smaller the uj/m ratio, 
the more likely to detect the particle closer to the horizon 
where the "effective gravitational potential" decreases its 
total energy. 



III. DETECTION OF MASSIVE PARTICLES 
OUTSIDE REISSNER-NORDSTROM BLACK 
HOLES 

Let us consider now the line element of a globally- 
hyperbolic spherically-symmetric static spacetime 



f{r)dt^ - h{r)dr^ - r^{d9^ + sin^ 9d(p^) , (3.1) 



where we assume that we are restricted to a region 
(/(r) > 0) where the Killing field {dt}^ is timelike. We 
may cast the positive-frequency solutions of the Klein- 
Gordon equation (□ + TO^)M2im(^) — the form 
{I eN, -I <m<l) 
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FIG. 4. The scattering potential is plotted for I = and 
I — 10. Full and dashed lines are related with Q = and 
Q — 0.9M, respectively. vanishes at the horizon and 

tends asymptotically to m^. We note that the larger the 
angular momentum, the higher the potential and the more 
difficult to find a particle with u/m < 1 far from the horizon. 
(We have assumed here mM = 2.) 



where i/j^iir) satisfies 



(3.2) 



with the following scattering potential 
fjjhd^fjh , l{l + l)f 



cff 



dr 



(3.3) 



(3.4) 



The label a was introduced to distinguish betwe en t he 
two independent solutions of the radial equation (|3.3|) . 

Outside a static black hole (r > r+ ) with mass M and 
charge Q (Q < M), we have /^^(r) = l/h^^{r), where 



/I^^(r) = (l-r+/r)(l-r_/r)>0 



(3.5) 



and r± = M ± a/M2 - Q^. (The label "RN" is a short- 
cut for Reissner-Nordstrom.) In this case, the scattering 
potential 



* cit 



RN 



2M 

r 



21 r 



2M 



2Q^ 

r.4 



1(1 + 1) 



(3.6) 



FIG. 5. The scattering potential is plotted for mM = 1/8 
and mM — 2. Full and dashed lines are related with Q = 
and Q = 0.9 M, respectively. We note that the larger mM, 
the more favorable to detect particles with uj/m < 1. (Here 
we have assumed / = 0.) 



vanishes at the horizon, r — r+, and tends to rn^ asymp- 
totically [see Fig. (^]. Thus, only outgoing" particles, 
a =—^, from the white hole horizon 7i~ will be able to 
comply Eq. (Eventually they will be totally re- 

flected to the black hole horizon Ti.^ .) Incoming particles 
from J'~ , a will necessarily have to > m. [The case 
Lu = m is somewhat subtler and will not be considered 
here (see Refs. [0 for the case lj — m ~ 0).] Let us 
note that the larger the I, the higher the centrifugal bar- 
rier of the potential [see Fig. (jj)]. As a conseque nce, par- 
ticles with / 3> 1 which are found to verify Eq. (LI) will 
be necessarily restricted to a close neighborhood around 
the horizon, where the gravitational field is strong enough 
to compensate the angular momentum effect. Hence, in 
order to find particles satisfying Eq. (LI) "relatively" far 
from the horizon, we will focus on particles with I = 0. 

For chargeless black holes, the local extrema of the 
scattering potential can be analytically determined. For 
the case I = and Mm < 1/4, we have 



1 



sm ■ 



cos ■ 



and r2 = — cos f , (3.7) 

TO 3 



where ri and r2 are associated with a local maximum 
and minimum (ri < r2), respectively, and ^ = tt — 
arctan[-y/l — 16(TOM)2/(4mM)]. For the case I = and 
Mm > 1/4 there are no local extrema [see Fig. (||)]. 
For Mm = 1/4, ri — r2 is an inflection point and 
V^\r=r^,r2 = 3toV4. For a black hole with M > 3Mq 
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and a particle with mass m > m^- , we have Mm > 10^^. 
This is why we will focus our simulations on values 
Mm > 1/4. An analytic investigation of the local ex- 
trema when Q ^ is not possible since dV^ /dr — 
results in a full fifth-order algebraic equation but it can 
be seen numerically [see, e.g., Figs. and (||)] that the 
presence of charge does not change significantly the form 
of the potential V^{r/r+). 

By using the coordinate transformation 



X = y + 



V- ln(?/ -y-)-y\ \n{y - y+) 

y- - y+ 



(3.8) 



where y = r/2M and y± = r±/2M, we can rewrite 
Eq. Q as 



'2„/,a 



dx"^ 



AMXTi4yir)WZi = ^M'lo'^Zi ■ (3.9) 



Thus, close to and far from the horizon, we can write the 
outgoing ip"i functions as > 0, i.e., Wmin — 0) 



j^-, -2tMux 

UJL 



(x < 0,|a;| > 1) 
(a;> 1) 

(3.10) 



and the incoming ones as (w > to, i.e., uj„ 



i^ujli^) ~ -^Lul i „-2iMCjx _|_ 'j^^^2iM{ux 



(x < 0,|a;| > 1) 
(a:» 1), 



where w 



m^ 



For iv>m, |7^;:;^|^|7^: 



(3.11) 
1^ and 



, l"?^!^ can be seen as reflection and transmission 
coefRcients, respectively, satisfying the probability con- 
servation equations 



+ - 1^.71' - 1 and \n::;f + ^ |Tjr|^ = i. 

LU to 



(3.12) 

For Lu < m, the outgoing modes fade exponentially far 
away from the horizon and = 1. One of the tests 

performed to guaranty the reliability of our codes was the 
verification that Eqs. ( 3.10 )-( 3.11 ) were satisfied (with 
the appropriate relations for 7?."; and 7^) along the nu- 
merical calculation. 

We normalize ip°^i such that u'^i^ are Klein-Gordon 
orthonormalized 



I I dT. n^' 

'St 



i" \7 11°' 



(3.13) 
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FIG. 6. We plot the detection rate F^^^qq as a function of 
the uio/m ratio for observers at different points ro. Full and 
dashed lines are related with Q = and Q — 0.9M, respec- 
tively. The damping in the detection rate is steeper for ob- 
servers far away from the horizon. (We have assumed here 
mM = 2.) 



'uj'l'm' 



= 0, 
(3.14) 

9AA' = ^ aa' ^ U'S mm' ~ ^^'^ '^^ the 

future-pointing unit vector normal to the volume element 
of the Cauchy surface S^. By using Eq. (3.9) to trans- 
form the left-hand side of Eq. (|t|) in a surface term, 
we obtain the normalization constants = (2uj)~^ and 



where here 6/ 



a: 



(2> 



^ (up to an arbitrary phase) . The incom- 



ing and outgoing modes are orthogonal to ea ch oth e r wit h 
respect to the Klein-Gordon inner product ( 3.13 )-( 3.14 ). 
This can be seen by choosing Ef — H^UJ^^ in Eq. ( 3.13 ) 
and recalling that ip^iix) and 'ip^iix) vanish on J^^ and 
Ti.~ , respectively. 

Next we expand a massive scalar field $(a;) in terms of 
positive and negative frequency modes: 



*(^) = E E E / '^^ Klm<lmi^) + H.C.] , 
a=^ 1=0 m=-l "''^min 

(3.15) 

where Wmin — when a and LOmm — ™ when a =^. 
As a consequence of orthonormalizing the normal modes 
with the Klein-Gordon inner-product, the annihilation 
'^Zim a^*^ creation dZj^ operators satisfy the simple com- 
mutation relations 
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FIG. 7. The probability density {dV^oo / dXd) / m is plotted 
for uj/ra — 0.4 and ui/m = 0.8. Full and dashed lines are 
related with (3 = and Q — 0.9Af , respectively. (Here we 
have assumed Mm = 2.) We note that the smaller the uj/m 
ratio, the closer to the horizon the particle should be. 



(3.16) 



for 



The Boulware vacuum |0) is defined by a^/jjjO) 
every a,u!,l and m. 

Let us assume an Unruh-DeWitt detector coupled to 
the massive scalar field as described by the interaction ac- 
tion (2^). The detection rate (detection probability per 
detector proper time) of particles with quantum numbers 
a, Lu, I and m (as defined by asymptotic fiducial observers) 
can be calculated (at the tree level): 



dim 



(3.17) 

where rip = F^{ll>) — const., we have chosen (3eq{^) = 
5{E — Eq) and we have tuned again the energy gap of 
the detector at each point in order to maximize the 
detection probability, namely, Eq—loI y/J^^jd). 

As in the case of the Rindler wedge, we may calcu- 
late the detection rate F"^;^ in the particular case where 
the massive state \auJolm) is defined by an experimen- 
talist lying on the precise location of the detector, i.e., 
1^0 = i^d- In principle, this would require that the scalar 
field were quantized with respect to a fiducial observer 
at To (in which case e~*"* in Eq. ( |3.2D would be replaced 
by e-'^^ with r = V/^^(ro)0- K is not difficult to 
see, howev er, that ^^^i^ can be obtained directly from 
Eq. ( 3.17 ) with the following replacements: rd — s- ro 



and UJ ^ f^^ifo) i^o- [One can check this str ategy 
in the Rin dler wedge by obtaining directly Eq. ( 2.12| ) 



from Eq. ( 2.11 ).] In Fig. (^ we plot F^qo a func- 



tion of ujo/m for observers at different locations. (We 
recall that in four-dimensional spacetimes cq is dimen- 
sionless, in contrast to the case where the spacetime is 
two-dimensional.) 

Now, let us define from Eq. (3.17) the (scalar) normal- 
ized probability density 



dvxJdVd ^ TZU^d)/ 1 rzu^')dv' 



(3.18) 



Here {dV"i^/dVd)dVd is the probabihty that a particle 
with quantum numbers (a, to, m) be found at some spa- 
tial point = (rd,9d,4'd) in a proper volume dVd = 
dXdr'^sm9dd9dd(j)d, where dXd = \/ h^^ {rd)drd. The 
proper radial distance can be integrated from the horizon 
r = r+ as a function of the radial coordinate: 



A(r) = / dr'^m^{r') 



r+ -M 



(3.19) 



Notice that \ — X{r) is a regular function except in the 
case of extreme Reissner-Nordstrom black holes [Q = 
M). From Eq. (3.18), we obtain the radial-distance prob- 
ability density 



dV"^, 



/in 



d\d /o+°°dA,'V?^^|^2,(r,')|^ 



(3.20) 



Note that here {dV"i^/dXd)d\d is the probabihty of de- 
tecting a particle with quantum numbers (a,a;,Z,m) in- 
side a she ll between and A^ -I- dXd- In Fig. (0), we use 
Eq. ( 3.20D to plot dV"ij^/dXd for particles with different 
uj/m ratios (and I — 0). We see that the smaller the oj/m 
the closer to the black hole horizon the particle will be 
in average. 

Next, in order to determine the mean particle radial 
distance, we calculate numerically [see Fig. (m] 



dXd XddVZim/dXd 



(3.21) 



Note that modes with Lo/m > 1 spreads over the whole 
space and thus (A^) diverges. This is qualitatively dif- 
ferent from the Rindler wedge case where the scattering 
potential grows unboundedly at the infinity what causes 
the normal modes to vanish asymptotically. As a con- 
sequence, at the Rindler wedge, (pd) is finite even for 
Lu/m > 1 [see Eq. (^T7|)]. 
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FIG. 8. (\d) is plotted for a —-^ and I — and com- 
pared with A™'"". We see, indeed, that (Ad) < A™'''' in the 
"high-frequency" regime uj > r^^ (i.e., at the right of the 
vertical broken line). (We have assumed here mM = 2 and 
Q = 0.) 



Now, by repeating our General Relativistic analysis 
performed in the Rindler wedge case in a more com- 
pact way, we obtain that according to (our fiducial) 
static asymptotic observers a particle with mass m at 
some point (r^, 9p, (f)p) with some four- velocity u'^ has to- 



tal energy lo = mu^{dt)^, = ™V/^^(?'p)\/l + I 
'^v/^^^V)- (The equality holds when the particle is at 
rest.) By inverting this relation, we obtain that a classi- 
cal particle with mass m and total energy lo (with lo < ni) 
will be at Vp < r^'^^, where 



M[l + ^1 - (1 - ujym^){Q/My 

1 — LO"^ jw? 



(3.22) 



(and Tp = r^'^^ only when the particle lies at rest). Be- 
cause the proper radial distance A = A(r) is a grow- 
ing monotonic function, Xp < A™'"'' = A(r™''''). In the 
"high-frequency" regime, u) ^ r^^, we expect as be- 
fore (Arf) < A™'*'' (see Sec. |l|). This is indeed verified 
in Fig. (|). 



IV. DETECTION OF MASSIVE PARTICLES AT 
RELATIVISTIC STARS 

Let us consider a static relativistic star with uniform 
density ao = const. The associated line element can be 
written as in Eq. (p.l[) (see Ref. B) with 



e(r - R) 



(4.1) 



and 



h'{r) = (1 - 2m(r)/r)-i 



1 - 2M— 1 e(i? - r) 



r I 



(4.2) 



where m(r) = 47r aor'^dr' . We are using here the label 
"s" to denote quantities associated with the star space- 
time. In these coordinates, r — R defines the star radius 
and M = (47r/3)aoi?'^ is the star total mass (which, of 
course, differs from the star total proper mass because of 
the binding energy contribution). Note that for sake of 
stability, R> Rc = 9M/4. 

The positive-frequency solutions can be cast in the 
same form given in Eq. ( |3.2[ ) where '(/'"; is defined by 
Eq. (|U) with /(r) and h{r) given in Eqs. (g^T]) and (^), 
respectively. The scattering potential can be obtained by 
using Eqs. (O) and (O) in Eq. (Q: 




-9F{R) 9 F{R) 1 
~ 2 V /i«(r) ~ h'ir) 



1 fl{l + l) 2 



1 



xe{R-r) + f'{r) 



f2M 1(1 + 1) , , , 

— + ^ .-^ ' +m^ ]eir-R) (4.3) 



where F{R) = 1 — 2M/R. At the star surface the poten- 
tial has a discontinuity: 



lim V,Ur) = FiR) 



1(1+ I) 
i?2 



M 
R^ 



In Fig. H we plot the scattering potential for different 
star parameters. (Note that the discontinuity of V^^ at 
r/R = 1 is very small for the parameters chosen in Fig. ^ 
to be self-evident.) It is interesting to observe that the 
gravitational field of the Earth is too faint to significantly 
infiuence the scattering potential for an electron-mass 
particle (m — 0.5 MeV). At the center where the influ- 
ence would be maximum, we would have for AI — G-lO^^'g 
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FIG. 9. The scattering potential is plotted for different 
R/Rc values. We have assumed here mM — 2. The denser 
the star, the larger the influence of the gravitational field. The 
first eigenfrequencies for these potentials are listed in Tab. |. 



10 '^^) which should 
^ obtained in the 



(and I = 0) that V^'g{r) = m^{l 
be compared with V^ff{r)\j^j^^ = m 
absence of gravity. Notwithstanding, for quasi-extreme 
stars {R w Rc) the scattering potential can be quite af- 
fected by the gravitational field (see Fig. ^) . Indeed it is 
possible to devise situations where the gravitational field 
of a star would be of some importance for particle-physics 
processes 

The normal modes with uj < m will be bounded 
states and Wmin = '-^i will be the lowest eigenvalue 
of the discrete frequency spectrum. Moreover, there 
will exist only one set of normalizable solutions and so 
we will suppress the "a" label hereafter. By defining 
dr = {2M)-^y^lHryjHr)dr, we can write Eq. dJ) as 



df'^ 



+ 4Af2l/4[f (r)]V^, = AMW^^i . (4.4) 



The function r{r) is a growing monotonic function. At 
the star core (and choosing properly the integration con- 
stant) we have f|r~o ~ (r/M)/(3^i^(i?) - 1) and the 
scattering potential becomes 



4Af2f2 



(4.5) 



where 
C 



2_G{Rf \{mMf , M3Z(; + 1) , Af^ Af3 



2M2 



R^G{R) R^G{R) i?3 

(4.6) 



FIG. 10. We plot the detection rate T^^qq for the eigenfre- 
quencies uo/m < 1 assuming observers at different points tq. 
We note that the damping in the detection rate is steeper for 
observers at weak gravitational potentials. (Here mM = 2 
and R/R^ = 1.2.) 



is a positive-definite constant and G{R) = — l + 3-\/F(i?). 
Thus, close to the star center, -0^^; is proportional 
to 2Mvfii{2Mvf) (since the solutions proportional to 
2Mvfni[2Mvf) are non- normalizable), where ji{x) is the 



spherical Bessel function and v 
far away from the star, V^^{r) s; 



\JlJ^ — . Moreover, 
Hence 



, B^i2Mvfji{2MvT) (r~«0) 
^ + e2'^^'= (f » 1) 



where we recall that = -vw^ — rn? . For uj > m, the 
modes are asymptotically free states and they can be 
normalized through Eqs. (3.13)-(3.14). Indeed we find 
^ujI = {2\/ujoj)~^ and \Cuji\ = 1. For iVi < m, the nor- 
mal modes fade exponentially and \Cc^-i\ = 0, where we 
are using here latin indeces to lable the discrete 

eigenfrequencies. The bounded states should be Klein- 



Gordon orthonormalized by imposing Eqs. (3.13)-(3.14) 
except for the fact that in this case 5aa' = SuiSnuii'^ujiLj.i ■ 
Thus, we find 



-f-oo 



^M2 



dr 1-01^,; (r) I 



4Mwf 



(4.8) 



The valid set of discrete eigenfrequencies Ui is determined 
numerically [|l0|. Shortly, the strategy consists in evolv- 
ing solutions with il^ujiiiO) = [see Eq. (4.7)] (and arbi- 
trary 'ipl,.i{f)\f=Q = const.) and search for the Wi's such 
that far enough from the star V'w,; ^ g-sAf^m^-i^f r 
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class 
^^min /'TI 


LJi/m 




uj-i/m 




mM = 2,R/Rc = 1.2 


0.264 


0.351 


0.456 


0.556 


0.648 


mM = 2,R/R^ = 2.0 


0.618 


0.679 


0.757 


0.826 


0.874 



TABLE I. We list the lowest eigenfrequencies \lo\, ...,bj^ 
for different star parameters. The lowest eigenfrequency uj\ 
is to be compared with the minimum classical energy i^Jjifn^. 



class 
'^min 
class A/f / d2 



regime i^—-' ^ M/R . By assuming that mAf 
M/i?^ = 0.07 m and M/R^ = 0.02 m for R/R^ 
R/Rc — 2, respectively. 



2, then 
1.2 and 



ruled by Eq. (4/7) for uji < m. In Tab. | we present 
the lowest eigenfrequencies {cji, 0^2, W3, 1^4} (normalized 
by m) for some star parameters. For sake of compari- 
son, we also show the minimum energy value djj^^jj*' for 
a classical particle obtained according to General Rel- 
ativity when it is static at the star center. Namely, 
<?f = my^^lr^Q = (m/2)(3Vl-2M/i?- 1). In 
the "high-frequency" regime, uj'^^'^^^ ^ M/R^, where the 
behavior of classical and quantum particles can be com- 
pared, we expect oji > ui^f^^ because of the extra intrinsic 
"kinetic energy" of quantum origin (see Tab. |) 



The massive scalar field can be cast as in Eq. (3.15) 



+00 ; 



$(a;) 



(a;) -HH.c.; 



(4.9) 



where the summation in a was suppressed and we in- 
tegrate over the free states [oj > m) and sum over the 
eigenfrequencies of the bounded states {uji < m). Here 



5aa' 



(4.10) 



where Saa' = (5;i'(5,nm'^(^ - t^') and 5aa' = 6w 5mm' ^<^cj' 
for free and bounded states, respectively. 

Next we couple [see Eq. ( |2.8| )] the massive scalar field to 
an Unruh-DeWitt detector characterized by the density 



of states (2.7). The amplitude (at the tree level) associ- 
ated with the detection of a bounded state with quantum 
numbers w^, / and m is A^fi^ = (0| (g) {E\Si\Eg) «> |wjm). 
The detection rate (detection probability per detector 
proper time) is, thus, 

r^ilm = -tot / dE PEoiE)\Atflm? 

Jo 

= 2cgny7^^^:^|l^„(e,,<^,)P, (4.11) 

^d 

where again the detector selectivity was chosen such that 
(i? I m (0)1 -Eg) = 1 and we have tuned the energy gap of 



0.4 



w/m=0.351 



^ ^cj/m = 0A56 
' ^ cj/m=0.556 
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FIG. 11. The probability density {dVujOO / d^d) /m is plot- 
ted for the three lowest eigenfrequencies when mM = 2 and 
R/Rc = 1.2 The smaller the u/m ratio, the closer to the star 
center the particle is. 



the detector as usuall y: E q = LUi /\/ f ^irg)- [Notice the 



resemblance of Eqs. ( 4.11 ) and ( 3.17 ), where we recall 
that dim(7/'^.i/'02;) = dim(\/a;).] (The dete ction of free 
states follows as previously [see, e.g., Eq. ( 2.10 )].) As 
in the Reissner-Nordstrom case, we can calculate the de- 
tection rate T^q;™ in the special case where the massive 
state \ujqItcl) is defined by an experimentalist lying on the 
detector, as sho wn in Fig. |l^. 

By using Eq. (4.11), we can define a normalized prob- 
ability density dV^jiim/ dVd analogously to Eq. (3.18), 
where dVd = d\ sin Oddcf) a and d\d = yjh''{rd)drd. 
The proper radial distance can be integrated from the 
center of the star as a function of the radial coordinate 
leading to 



^ ^ arcsin[V2MrVit;3] ^^^ _ 



arcsin[^2M/i?] 
^2M/i?3 



+Mln 



r + r^Jl - 2M/r - M 
R + R^l - 2M/R - M 

<d[r-R] . 



r\ l- 



2M 




(4.12) 



The radial-distance probability density dVu,iim/dXd 
then calculated analogously to Eq. ( 3.2C| ). In Fig. (|l]^ 
we plot dVi^-im/dXd for particles with different uji/m ra- 
tios. The mean particle radial distance (Ad) is defined 
analogously to Eq . (|3.21 ) and was calculated numerically 
as shown in Fig. ([l^ 
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FIG. 12. (Xd) is plotted for various eigenfrequencies as- 
suming I = 0, niM = 2 and R/R^ = 1.2. We note 
that (Ad) < A™™, as expected, in the "high-frequency" 
regime a;'='''=7m > {M/R^)/m. (Here cu^fn/m = 0.264, 
oJi/m = 0.351 and (M/R^)/m = 0.07.) 



Now, by taking a procedure analogous to the one used 
in the black hole case [see Eq. ( 3.22| )], we find that a 
classical particle with total energy to < m (as measured 
by static asymptotic observers) will be at rp < r^'^^, 
where 




m 



(4.13) 



Because the proper radial distance A = A(r) is again a 
growing monotonic function [see Eq. ( 1.12| )] , we have once 
more Ap < \^^^ = X{r^^^) . This is compared with (A^) 
in Fig. (H). We see that (Ad) < X^p'^^, as expected, in the 
"high-frequency" regime, oj'^^'^^ ^ M/E?. 



V. CONCLUSIONS 

The detection rate for massive particles with total en- 
ergy E < m(? was calculated for the Rindler wedge, black 
holes and star spacetimes. The mean particle positions 
were calculated and shown to be in qualitative agree- 
ment with General Relativity predictions in the "high- 
frequency" regime. The observation process is natu- 
rally taken into account in the formalism by consider- 
ing Unruh-DeWitt-like detectors. In this way, we clearly 
define what we mean by "observing a particle" . 



Although it is possible, in principle, to measure parti- 
cles satisfying Eq. (1.1) even at Earth, these events are 
so rare in such weak gravitational fields that do not have 
to be considered for practical purposes. This can be seen 
from the analysis of the scattering potential [see discus- 
sion below Eq. (1.3)]. 

Notwithstanding the consideration of massive particles 
with arbitrary small total energy is just fundamental for 
a comprehensive understanding of some phenomena oc- 
curring in highly curved spacetimes (e.g., black holes) 
and accelerated frames. For example, the weak decay of 
uniformly accelerated protons in the Minkowski vacuum 
which is described by the p ^ n v channel (at the 
tree level) in an inertial frame was shown to be repre- 
sentable by the combination of the following three chan- 
nels: p e~ -^nv^pp—^ne^ and p e~ 9 n , 
in the proton coaccelerated frame. According to this de- 
scription, the proton would decay by the absorption of 
a Rindler e~ and/or D from the FuUing-Davies-Unruh 
thermal bath "attached" to the proton proper frame. 
However, in order that inertial and noninertial frame 
descriptions lead to the same observables (e.g. proper 
lifetime) , we emphasize that the energy spectrum of the 
(massive) Rindler e~'s and e+'s must he cu G [0, -l-oo) 
(despite Rindler e^'s and e+'s have mass m « 0.5 MeV, 
as usually) Perhaps an even more extreme case con- 
cerns the importance of zero-energy particles to under- 
stand the radiation emitted from uniformly accelerated 
charges according to comoving observers. Namely, the 
emission of a (usual) Minkowski photon with transverse 
momentum from a uniformly accelerated charge as 
described by inertial observers corresponds, according to 
coaccelerated observers, to either the absorption from or 
the emission to the Fulling-Davies-Unruh thermal bath of 
a zero-energy photon with the same transverse momen- 
tum kj_. For obvious reasons, in Ref. [|ll| the photons 
considered were supposed to be massless but the same 
conclusion in terms of zero-energy particles would hold 
even if we considered that the photons were massive. 
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